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STATISTICAL  MECHANICS  OF  A 
CONTINUOUS  MEDIUM  (VIBRATING 
STRING  WITH  FIXED  ENDS) 

J.  KAMPfi  DE  FfiRIET 

UNIVERSITE  DE  LILLE  AND  UNIVERSITY  OF  MARYLAND 


1.  For  the  past  twenty  years  there  has  been  much  discussion  in  fluid  mechanics 
about  the  statistical  theory  of  turbulence.  It  was  J.  Boussinesq  (1872)  and  O. 
Reynolds  (1895),  in  their  pioneering  work,  who  expressed  the  idea  that  the  turbu¬ 
lent  velocity  fluctuations  of  a  fluid  were  much  too  complicated  (changing  rapidly 
from  one  time  and  one  point  to  another)  to  be  known  in  all  their  details;  we  must 
be  satisfied  to  study  only  some  convenient  mean  values.  The  systematic  use  of 
statistical  methods  has  led,  since  1930,  to  very  important  results,  in  the  fundamen¬ 
tal  works  of  Sir  Geoffrey  Taylor,  Th.  von  Karman  and  A.  N.  Kolmogoroff,  for 
instance. 

But,  if  we  look  carefully  at  all  the  results  so  far  obtained,  we  see  at  first  glance 
that  they  are  not  at  all  in  the  same  close  relation  with  the  theoretical  equations  of 
fluid  mechanics  as  the  classical  statistical  mechanics  bears  to  the  Hamilton-Jacobi 
equations  for  a  dynamical  system  having  a  finite  number  of  degrees  of  freedom. 
For  such  a  system  the  main  features  of  statistical  mechanics  are  the  following: 

(a)  definition  of  the  phase  space  12;  every  state  of  the  system  is  characterized 
by  a  point  co  £  12,  whose  2k  coordinates  are  the  Lagrangian  parameters  qi,  .  .  .  , 
and  the  conjugate  moments  Pi,  ,  pk. 

(b)  proof  of  a  uniqueness  theorem:  starting,  at  the  initial  time  t  =  0,  from  a  given 
initial  state  o>,  all  the  subsequent  or  prior  states  Tt  co,  where  t  varies  from  —  oo  to 
+  oo ,  are  perfectly  well  determined  and  describe,  in  12,  a  curve,  the  trajectory  rw. 

(c)  definition  of  a  measure  in  12,  invariant  under  the  transformation  Tt  co.  (Liou- 
ville’s  theorem.) 

(d)  proof  of  the  ergodic  theorem  (existence  of  time  averages  computed  along  a 
trajectory  rw);  the  time  averages  and  statistical  averages  can  be  validly  con¬ 
sidered  as  equal  only  in  the  special  case  in  which  the  transformation  Tt  co  of  12  into 
itself  is  metrically  transitive. 

As  we  have  pointed  out  in  [4],  at  the  time  being,  all  the  work  done  in  what  is 
called  statistical  theory  of  turbulence  deals  only  with  some  kind  of  time  averages; 
so  far  as  we  know,  a  convenient  phase  space  S2  has  never  been  introduced,  in  order 
to  represent  all  the  possible  states  co  of  a  fluid;  this  phase  space  must  be,  as  we  have 
shown  in  [5],  a  function  space,  quite  easy  to  define.  But,  so  little  is  known  about 
the  solutions  of  the  nonlinear  equations  of  fluid  mechanics  (for  a  perfect  as  well  as 
for  a  viscous  fluid),  that  immediately  after  the  definition  of  S2,  we  are  stopped  by 
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the  lack  of  a  uniqueness  theorem  and,  especially,  by  the  lack  of  any  method  leading 
to  the  definition  of  an  invariant  measure  in  12. 

To  point  the  way  to  what  has  to  be  done  in  order  to  build  a  real  statistical 
mechanics  of  turbulence,  it  has  seemed  worth  while  to  us  to  sketch  a  statisti¬ 
cal  mechanics  for  a  continuous  medium,  much  simpler  than  a  fluid.  For  a  vibrating 
string,  due  to  the  linear  character  of  the  equations,  the  problem  can  be  completely 
solved;  in  [5]  we  have  treated  the  case  of  the  infinite  string;  we  want  to  discuss  here 
the  string,  with  fixed  ends,  which  has  much  more  interesting  features  than  the  in¬ 
finite  one. 

2.  In  the  manner  classical  since  Bernoulli,  Euler  and  D’Alembert,  let  us  consider 
a  vibrating  string  as  a  material  curve  moving  in  a  plane  0 xy;  when  at  rest  the  string 
is  the  straight  line 

0  £x£l,  y  =  0. 

Every  point  x  of  the  string  moves  parallel  to  the  Oy  axis  and  we  represent  by 
y(x,  t )  its  displacement  at  the  time  /. 

We  will  make  the  four  following  assumptions: 

(2.1)  y  (a,  t)  is  a  continuous  function  of  (»,  t)  , 

(2.2)  y  (0,  l)  =  y  (l,  t)  =  0  for  every  t , 

(2.3)  y  (x,  l)  has  continuous  derivatives  up  to  the  second  order  , 

(2.4)  y  (x,  t )  is  a  solution  of  y'J  ( x ,  t)  =  yy  (#,  t)  . 

For  the  sake  of  brevity,  let  us  put 

(2.5)  p  (x,  t)  =  y'  (x,  t) , 

(2.6)  v  (x,  t)  =  y't  (x,  t)  , 

(2.7)  j  (x,  t)  =  y,"  (x,  t)  ; 

v  and  j  represent  respectively  the  velocity  and  the  acceleration  of  the  point  x  of  the 
string  at  the  time  t. 

The  motion  of  the  string,  from  t  =  —  °°  to  t  =  +  <» ,  is  completely  determined 
by  the  knowledge  of  the  displacement  and  the  velocity  for  every  point  x  at  any 
given  time  t. 

In  fact  let  us  introduce  two  functions  of  x,  y(x)  and  i>(^),  such  that: 

(2.8)  the  displacement  y  ( x )  is  a  continuous  function  for  0  ^  x  ^  l , 

(2.9)  y  (0)  =  y  (/)  =  0  , 

(2.10)  the  first  two  derivatives  y'(x)  and  y"(x)  exist  and  are  continuous  ; 

(2.11)  the  velocity  v  (#)  is  a  continuous  function  for  0  ^  x  ^  l , 

(2.12)  v  (0)  =  v(l)  =  0, 

(2.13)  the  first  derivative  v'  ( x )  exists  and  is  continuous  . 
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From  these  functions,  we  derive  a  function/^)  defined  for  —  »  <  x  <  +  °°  by : 


(2.14) 

nx)=y'(x)  +  vMt 

O^x^l, 

(2.15) 

/(*)=?'(-*> -■(-*) , 

-l^x^O, 

(2.16) 

f(x  +  2l) -/(*). 

It  is  obvious  that/(#)  is  a  continuous  function  of  x,  having  a  continuous  deriva¬ 
tive  f'(x),  and  that 

(2.17)  P  f(x)dx  =  0. 

The  unique  displacement  function  y(x,  t),  satisfying  (2.1)  to  (2.4)  and  such  that 
y  (x,  t)  =  y(x)  ,  v  (x,  t)  =  v  (*) 
is  given  by  the  classical  formula 

(2.18)  y{x,t)=  f  f(l-T+s)ds, 
from  which  we  obtain  immediately 

(2.19)  p  (x,  t)  =  f  (t  —  T  +  x)  +/(/  -  T  —  x) , 

(2.20)  v  (x,  t)  =  /  (/  —  r  +  x)  —  f  (t  —  t  —  x)  , 

(2.21)  j  (x,  t)  =  f  (t  —  t  +  x)  —  f  (t  -  t  -  x) . 

3.  According  to  the  preceding  well  known  results,  we  shall  take  as  phase  space  12 
the  following  function  space:  12  is  the  set  of  all  real  valued  functions  f(s)  of  a  real 
variable  s,  —  °°  satisfying  the  four  conditions: 

(3. 1 )  x  /  (s)  is  periodic,  /  (s  +  21)  =  f  (s)  , 

(3.2)  f  (s)  is  continuous  , 

(3.3)  /  (s)  has  a  continuous  derivative/'  (s)  , 

(3.4)  f  /(  s)  ds  =  0  . 

Every  function  f(s)  satisfying  these  assumptions  is  called  a  point  w  of  the 
space  12. 

Every  state  of  the  vibrating  string,  being  well  defined  by  the  knowledge  of  the 
displacement  function  y(x)  and  the  velocity  function  v(x),  is  represented  in  the 
phase  space  12  by  the  point  w  corresponding  to  the  function  f(s)  computed  accord¬ 
ing  to  (2.14),  (2.15),  and  (2.16).  Conversely,  from  every  point  w  £  12  we  can  de¬ 
duce  a  displacement  function  y(x)  and  a  velocity  function  v(x)  by  solving  (2.14) 
and  (2.15): 

(3.5)  y  (*)  =  f  f  (s)  ds  , 

J  —  X 
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hence  to  every  co  £  12  corresponds  a  completely  determined  state  of  the  string. 

There  is  thus  a  one  to  one  correspondence  between  the  states  of  the  string  and 
the  points  of  the  phase  space  12. 

We  may  note  that  12  is  obviously  a  separable  Banach  space,  if  we  define  the 
norm  by:1 

(3.7)  M|  =  max  |  /($)  | . 

0=8=21 

4.  Let  us  now  introduce  a  one  parameter  Abelian  group  of  one  to  one  transforma¬ 
tions  of  12  into  itself : 

(4.1)  CO~^TtCO,  —co  <  /  <  -j-  co  ; 

<jj  corresponding  to  f(s),  we  define  its  transform  Tt co  as  the  point  corresponding 
to  /(/  +  s). 

It  is  clear  that  for  every  co  £  12, 

(4.2)  TV  (7>«)  =  Tr  ( 7>co )  =  7Y+<"co , 

(4.3)  70co  =  co  , 

(4.4)  (r<)~1co=  T_t  co, 

(4.5)  T t+2ico  =  7)  co  . 

With  the  definition  (3.7)  of  the  norm,  the  transformation  is  isometric: 

(4.6)  M-IMI. 

If  we  consider  a  fixed  point  co  £  12,  the  set  of  its  transforms  71  co,  when  t  varies 
from  —  co  to  +  oo ,  is  a  closed  curve,2  which  we  shall  call  the  trajectory  rw.  The 
trajectory  T0  corresponding  to  the  point  co  =  0  contains  only  that  point;  except  for* 
this  particular  trajectory  T0,  we  can  establish  a  one  to  one  correspondence  between 
the  points  of  the  trajectory  r„  and  the  points  of  a  circle  C  of  radius  //-tt;  let  us  fix  the 
position  of  a  point  on  C  by  the  length  y  of  the  arc  ending  at  this  point, 

(4.7)  Ogy<2  Z; 
then  we  obtain  the  one  to  one  correspondence 

7— >  TyCO  . 

Thus  it  is  clear  that,  when  t  varies  from  —  oo  to  +  °° ,  the  transforms  Tt  co  of  co  de¬ 
scribe  the  trajectory  infinitely  many  times. 

1  Usually  for  the  Banach  space  of  the  continuous  functions/^)  having  a  continuous  derivative 
f'(s),  the  norm  is  taken  as 

IMI  =  max  \f(s)  |  -1-  max  |  f(s)  |  , 

0£s<2l  Q£s<2l 

but  here,  because  of  (3.4),  the  second  term  need  not  be  introduced. 

2  The  name  “curve”  is  justified  by  the  fact  that  the  continuity  condition 

lim  ||T t+h  «  -  7\«H  =0 
0 


is  satisfied  for  every  t. 
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Using  the  geometrical  language  we  have  introduced,  we  can  restate  the  unique¬ 
ness  theorem  of  section  2,  in  the  following  terms:  when  t  varies  from  —  <»  to  + 00 , 
all  the  states  of  the  vibrating  string,  corresponding  to  given  initial  conditions,  are 
represented  in  the  phase  space  ft  by  the  unique  trajectory  T*,  which  passes  through 
the  particular  point  co  representing  the  initial  state. 

The  fact  that  Ta  is  a  closed  curve  means  that  after  a  period  21  the  state  of  the 
string  is  again  the  initial  state:  the  motion  is  periodic. 

The  singular  trajectory  T0  corresponds  to  the  case  in  which 

y  (x)  =  0 ,  r  (*)  =  0; 

the  string  is  always  at  rest  and  the  unique  point  u  =  0  represents  all  its  states. 

6.  In  what  follows  we  shall  represent  by  a  symbol  like 

F  [o>  |  Xl,  .  .  .  ,  Xn] 

a  functional  defined  on  ft,  depending  on  n  parameters:  to  every  point  co  £  ft  and  to 
every  set  of  values  of  the  parameters  Xi,  .  .  .  ,  Xn  corresponds  a  real  number. 

Of  special  interest  are  the  linear  functionals;  in  our  particular  space  ft,  they  are 
given  by  the  F.  Riesz  theorem: 

(5.1)  L  [w|  Xi,  .  .  .  ,XJ  =  f  /($)<*$($,  Xi,  ...  ,Xn), 

‘/0 


where  <i>(s,  Xi,  .  .  .  ,  X„)  is  supposed  to  be  a  function  of  bounded  variation  of  5  for 
every  set  of  values  of  the  parameters  Xi,  .  .  .  ,  Xn. 

We  shall  have,  in  particular,  to  consider  the  linear  functionals: 


(5.2) 

F[«|»]  =  f’/(,s)ds, 

0  gxgl. 

(5.3) 

=/(*)  +/(-  *)  > 

(5.4) 

V[a\x]  =f(x)-f(-x). 

Obviously,  according  to  (2.18),  (2.19),  and  (2.20),  along  the  trajectory  ru 
(where  w  represents  the  state  of  the  string  at  t  =  0) : 

(5.5)  Y[Tt<a\x]  =  y  (x,  t) , 

(5.6)  P[Ttu\x]  =  p  (x,  t) , 

(5.7)  V  [rtw|a;]  =  v  (*,  t) . 

6.  Now  we  have  to  define  an  invariant  measure  m  in  the  phase  space  ft;  an  in¬ 
variant  measure  is  a  set  function  m(E)  defined  for  every  subset  £  of  ft  belonging  to 
a  <r-subalgebra  S  of  the  Boolean  algebra  of  all  subsets  of  ft  (Borel  field),  with  the 
following  properties: 

(6.1)  0  m  ( E )  ^  1  for  every  £  £  S  ; 

(6.2)  m  (ft)  =  1  ; 

(6.3) 

L  n  J  n  / 
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for  every  finite  or  countable  sequence  of  disjoint  sets  En  €  S,  Ej  n  Ek  =  0  for 
j  k; 

(6.4)  if  E  6  S  ,  then  TtE  6  S  for  every  t  and  m  ( TtE )  =  m  (E)  . 

We  can  obtain  the  measure  m  by  a  well  known  process  due  to  Kolmogoroff  [6]. 
We  start  with  cylinder  sets;  a  cylinder  set  in  0  is  the  set  of  all  <*>  corresponding 
to  functions /($)  satisfying  a  condition  such  as 

(^■5)  /  ($0  <  £1,  •  •  •  ,/  CO  <  £«  ; 

we  may  take  for  n  all  the  integer  values  1,  2,  .  .  .  ,  for  ($1,  .  .  .  ,  sn)  any  set  of  n 
different  real  numbers  and  for  (£1,  ...,£„)  any  set  of  real  numbers. 

Now  we  define  the  measure  of  the  cylinder  sets,  by  choosing  an  arbitrary  se¬ 
quence  of  distribution  functions3 

(£l)  >  (£l,  £2j  $2  ^l)  j  •  •  •  y  (£l>  •  •  •  >  £n>  ^2  ^1.  •••»£»  -^l)  >  •  •  •  > 

and  by  putting4 

(6.6)  [/  (si)  <  £i (s„)  <  £„] 

CM 

(£l.  •  •  *  >  (nj  *^2  -^1)  •  •  •  y  S n  ^l)  • 

This  being  done,  we  have  only  to  extend  the  definition  of  the  measure  to  the 
er-subalgebra  5  generated  by  the  cylinder  sets;  it  is  known  that  this  extension  is 
unique  [1],  [2]. 

The  invariance  of  the  measure  m  thus  defined  is  a  consequence  of  the  fact  that 
the  distribution  functions  do  not  depend  on  the  Si,  .  .  .  ,  sn,  but  only  on  the 
differences  s2  —  Si,  .  .  .  ,  sn  —  Si,  which  are  invariant  under  any  translation,  thus 

HE  If  (*  +  <  £i  >•••>/  (<  +  sn)  <  £n]  =  /  (sx)  <  £i ,...,/  ( sn )  <  £n]  ; 

60  CO 

now  it  is  quite  evident  that 

(Ji)  <  $1}  •  •  •  t  /  CO  <  £n]  =  +  Sl)  <  £l>  •  -  -  f  (t  +  $n)  <  £n]  , 

O)  u 

which  establishes  the  invariance  of  the  measure  of  the  cylinder  sets  under  the  trans¬ 
formation  Tt  u. 

*  As  is  well  known,  these  functions  must  satisfy  a  set  of  consistency  conditions: 
dpi(£i)  is  a  nondecreasing  function  of  &,  such  that 

42i(-c°)  =  o,  <£i(4®)  -  i ; 

£2,  52  —  Si)(s2  —  5i  0)  is  a  nondecreasing  function  of  (&,  £2)  such  that 
4**(— fa.  52  —  Jl)  =  —°°,S2  —  Si)  =  0, 

^2(4  00  >  $2,  52  ~  5i)  = 

^2(fl,  4  00  ,  52  —  Si)  =  ^i(£l), 

4^(+  °°>  4  ®,  52  -  5i)  =  1; 
and  so  on  (for  the  details  see  [6]). 

4  By  the  symbol  E[P]  we  mean  the  subset  of  12,  for  every  point  of  which  a  given  condition  P  is 

W 

'  satisfied. 
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But,  with  this  method5  the  difficulties  begin  when  we  try  to  take  account  of  the 
fact  that  12  is  not  the  space  of  all  real  functions,  but  only  of  functions  subject  to 
the  conditions  (3.1)  to  (3.4);  these  conditions  impose  on  the  distribution  func¬ 
tions  4^2,  •  •  •  some  particular  properties;  for  instance,  as  we  shall  see  later, 
in  order  to  satisfy  (3.4)  we  must  have 

(6.7)  =0. 

•/ — CO 

But  it  does  not  seem  at  all  easy  to  find  directly  all  the  conditions  resulting  for 
4?i,  •  •  •  ,  Qn,  •  ■  •  from  (3.1)  to  (3.4). 

4  To  avoid  this  difficulty,  we  shall  take  another  course,  making  sure  from  the  be¬ 
ginning  that  every  function  corresponding  to  a  point  co  of  12  does  in  fact  satisfy 
(3.1)  to  (3.4);  this  is  quite  easy  using  the  representation,  due  to  Norbert  Wiener, 
of  a  random  function,  for  the  particular  case  of  Brownian  motion. 

Let  us  take  an  abstract  measure  space  A ,  for  which  the  measure  n  is  such  that 

M  (it)  -  1  • 

Let  us  consider  a  real  valued  function  f(s,  a)  defined  over  the  product  space 
R  X  A,  where  R  means  the  set  of  all  real  numbers  —  ®  <  s  <  +<».  We  sup¬ 
pose  that  for  every  a  6  A,  as  a  function  of  s,  f(s,  a)  belongs  to  some  given  func¬ 
tion  space  12  (for  instance,  the  space  of  all  bounded  functions,  or  of  all  continuous 
functions,  etc.);  further  we  suppose  that  for  every  real  s,  as  a  function  of  a, 
f(s,  a)  is  measurable  (with  respect  to  the  measure  /z). 

When  a  describes  A,  the  set  of  all  functions  f(s,  a)  is,  in  general,  only  a  subset 
12*  of  the  function  space  12;  for  every  E*  c  12*  we  shall  denote  by  A(E*)  the  set 
of  all  a  whose  images  are  in  E*  (it  must  be  noted  that  the  correspondence  between 
A  and  12*  is  not  necessarily  one  to  one). 

For  any  subset  E  c  ft  we  shall  say  that  it  is  measurable  if  A[E  n  12*]  is  meas¬ 
urable;  if  this  is  the  case,  we  take  for  the  measure  of  E,  by  definition, 

(6.8)  m  (E)  =  fx  [A  (E  n  12*)] . 

It  can  be  proved  [3]  that  the  set  function  thus  defined  is  effectively  a  measure  on 
a  subalgebra  5  of  subsets  of  12. 

From  (6.8)  we  get 

(6.9)  m  (12)  =  m  (12*)  =  M  [A  (12*)  ]  =  1  ; 
for  any  E  such  that  12*  c  E 

m{E)  =  1  , 

and  for  any  E  such  that  12*  n  E  =  0 

m(E)  =  0 . 

8  Every  step  could  be  interpreted  in  the  language  of  probability  if  we  assume  that  we  pick  at 
random  a  point  w  in  Q,  the  probability  being  defined  by 

Pr[u£E\  =  m(E)  ,  E  6  5  ; 

the  definition  of  the  measure  that  we  have  sketched  corresponds  exactly  to  the  definition  of  a  ran¬ 
dom  stationary  function  of  s. 
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As  a  byproduct  of  this  definition  of  m,  we  obtain  the  distribution  functions,  by 
observing  that  the  set 

£  =  E  [/  (*l)  <  £1  ,•••>/  (*»)  <  £«] 

03 

has  the  same  measure  as  the  subset  of  £2* 

E  n  £2*  =Ef/0i,  a)  <  £1,  .  .  .  ,/(*»,  a)  <  , 

CO 

and  the  measure  of  this  set  is  simply  the  measure  of  the  subset  of  A  for  which 
/  ($1,  a)  <  £1,  .  .  .  ,  /  0„,  a)  <  £n  ; 

but  we  need  no  longer  worry  about  special  conditions  for  the  distribution  func¬ 
tions,  being  sure  a  priori  that  all  the  functions  f(s,  a)  we  have  introduced  possess 
just  the  desired  properties. 

In  order  to  define  an  invariant  measure  m  in  the  phase  space  £2  of  the  vibrating 
string,  it  is  very  convenient  to  take  for  A  a  product  space 

A  =  BXC, 

where  C  is  the  circle  already  considered  in  (4.7)  and  B  is  any  measure  space;  we 
take  on  C  the  ordinary  Lebesgue  measure  and  we  call  v  the  measure  defined  on  B; 
the  measure  n  on  A  is  then  defined  in  the  ordinary  way  for  a  product  space. 

Now  let  us  define  on  A  a  real  valued  function 

g  (ft  7) ,  0£B,  y^C, 


with  the  following  properties: 

(6.10)  for  every  7,  g  is  a  measurable  function  of  0  ; 

(6.11)  for  every  /3,  g  is  periodic  in  7,  g  (0,  7  +  21)  =  g  (0,  7)  ; 

(6.12)  for  every  0,  g  is  a  continuous  function  of  7  ; 

(6.13)  for  every  0,  g  has  a  continuous  derivative  gy  ( 0 ,  7)  ; 

(6.14)  for  every  0, 

f  g(0,  7)^7  =  0. 

•/o 


If  we  put 
(6.15) 


f(s)  =  g  ( 0 ,  s+  7)  , 


0£  B,  7 e  c, 


the  set  £2*  of  all  these  functions  f(s)  is  obviously  a  subset  of  the  phase  space  £2, 
since  the  conditions  (6.11)  to  (6.14)  imply  the  conditions  (3.1)  to  (3.4). 

The  meaning  of  0  and  7  is  clear;  the  set  of  points  co,  corresponding  to  a  given 
value  of  0,  and  to  all  the  values  of  7  in  the  interval  0  ^  7  <  22is  simply  a  trajectory 
T;  thus  £2*  is  really  a  set  of  trajectories,  each  value  of  0  defining  a  trajectory.  It 
follows  immediately  that  £2*  is  transformed  into  itself  by  Tt: 

(6.16)  Tt £2*  =  £2*  . 
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By  the  method  just  described  above,  starting  from  the  space  A  and  the  functions 
g(@,  s  +  7)  we  define  a  measure  m  in  12  by  the  measure  n  in  A ;  let  us  observe  that 
the  transformation  from  to  Tt «  corresponds  to  the  change  of  g(fi,  5  +  7)  into 
g(P,  t  +  s  +  7) ;  this  means  that  the  transformation  Tt  of  £2  into  itself  induces 
in  A  a  very  simple  transformation:  a  rotation  of  amplitude  xt/l  of  the  circle  C. 
The  Lebesgue  measure  on  C  being  invariant  under  a  rotation,  it  follows  that  the 
measure  n  on  A  is  invariant,  and  further  that  the  measure  m  on  ft  is  also  invariant. 
Thus  the  measure  m  fulfills  all  the  necessary  conditions  (6.1)  to  (6.4). 

The  invariance  of  the  measure  m  under  the  transformation  Ttw  of  the  phase 
space  ft  into  itself  is  just  the  equivalent  of  Liouville’s  theorem  in  classical  statistical 
mechanics;  it  is  a  necessary  condition  for  the  introduction  of  probability.  The  full 
sequence  of  all  the  states  of  the  vibrating  string  is  completely  determined,  accord¬ 
ing  to  the  uniqueness  theorem,  by  its  state  at  any  given  time  t.  A  definition  of  the 
probability  of  picking  at  random  a  state  of  the  string  makes  sense  only  if  it  gives 
the  same  value  for  the  probability  along  a  trajectory  T,  independently  of  the  time  t 
of  our  choice.  We  can  now  make  a  perfectly  coherent  definition  of  the  probability 
of  picking  at  random  a  state  of  the  vibrating  string  by  setting 

(6.17)  Pr  [co  6  E\  =  m  ( E )  for  every  E  £  S  , 

because  along  the  trajectory  r„  we  have  at  any  time  t 

Pr  [w(  TtE]  =  m  (£)  . 

7.  We  shall  say  that  a  functional  F[w],  defined  over  the  phase  space  ft,  belongs  to 
the  class  L,  if  its  Lebesgue  integral,  computed  with  the  measure  m,  exists. 

For  every  functional  F[«]  6  L,  we  define  the  statistical  average  by 

(7.1)  fR=  f  F  [co]  dm  =  fj?[u]dm. 

*/G  Q 

If  we  observe  that,  on  ft*,  any  functional  could  be  put  in  the  form 


(7.2) 

F  [«]  =  4>  (0,  7)  , 

we  get 

F  [co]  =  (p,  7)  dn  ; 

(7.3) 

since  A  is  a  product  space,  we  have  by  the  Fubini  theorem 

(7.4)  F[w]  =  ®  ^  7>  d'y\dv' 

Noting  that 

(7.5)  F[Ttco \  =  4>03,/+  7) 
and  that 

(7.6)  f  4>(/3,  *  +  7)  dy  =  f  4>  (0,  7)  dy  , 

*/o  ^0 

we  get  the  following  result: 

(7.7)  F[Ttu}  =  F  [col  for  any  given  t  . 
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As  in  classical  statistical  mechanics,  we  have  also  to  consider  a  second,  quite  dif¬ 
ferent,  type  of  average,  the  time  average;  this  average  is  equal  to  the  mean  of  all 
the  values  F[Tt  <o]  along  a  given  trajectory  ru: 

(7.8)  M.FM-*  Hm  -1=  fTF[T,»]dl. 

T — >  -f-  00  /l  •/ — 

From  (7.5)  we  get 

MtF  [«]  =  lim  f  4>(|8,  t  +  y)  dt\ 

T — >+00  LI  J — T 

thus,  for  every  functional  F[co]  £  L  and  for  every  trajectory  T„,  the  time  average 
exists  and  has  the  value 

(7.9)  Jf.FM  =-tjf!W  I)dy, 

this  result  is  the  particular  form  of  the  ergodic  theorem  for  the  statistical  mechanics 
of  the  vibrating  string. 

Any  time  average  Mt  F[<o]  is  itself  a  functional,  say  M[u>\ ;  this  functional,  being 
independent  of  7  according  to  (7.9),  is  invariant  along  every  trajectory  T: 

(7.10)  M[Tta]  =  M[ co]  . 

Combining  (7.4)  and  (7.9),  we  get 

(7.11)  FM  =  fBM.F\»]dv, 

proving  that,  except  perhaps  for  some  particular  functionals  (those  for  which 
MtFM  =  const,  for  every  «)  or  for  some  particular  choice  of  the  measure  m 
the  two  kinds  of  averages  are  not  equal;  this  result  means  that  in  general  the  trans¬ 
formation  Tt  of  ft  into  itself  is  not  metrically  transitive;  there  are  subsets  E  of  S2, 
with  a  measure  different  from  0  or  1,  invariant  under  the  transformation  Tt. 

We  obtain  two  interesting  examples,  by  taking  for  the  measure  space  B  the  two 
typical  cases: 

(a)  B  is  a  finite  or  a  countably  infinite  set  of  points  /3„;  we  define  a  measure 
v  in  B,  by  assigning  to  every  point  /3„  a  weight  vn  ^  0,  in  such  a  way  that 

y  vn  =  1  . 

n 

Then  the  measure  of  any  subset  Bi  of  B  is  simply  the  sum  of  the  vn  corresponding 
to  the  /3„  €  B\. 

In  this  case 

^[m]  ?)  <*7, 

and  for  the  trajectory  corresponding  to  /5„ 
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(b)  B  is  the  interval  0  ^  /3  ^  1  and  v  is  the  ordinary  Lebesgue  measure;  then 

7)  dyjdff, 

MtFM  =Ylf„ 

It  is  obvious  from  these  two  examples  that,  in  the  statistical  mechanics  of  the  vi¬ 
brating  string  with  fixed  ends,  it  is  impossible  in  general  to  replace  a  time  average 
by  a  statistical  average.  Yet  there  is  one  remarkable  exception:  for  any  linear 
functional  L[ w]  the  statistical  average  and  the  time  average  are  both  equal  to  zero: 

(7.12)  L  M  =  MtL  M  =  0  . 

According  to  (7.11)  it  is  enough  to  prove  that  >f/(L[a>]  =  0;  here  the  formula 

(5.1)  gives  for  any  linear  functional 

L[  co]  =  f  g(l S,  s  +  y)  d<l>(s)  , 

*/o 

whence 

*‘LM=7lfU‘  g(P>  s  +  7)  d<j>(s)  Jd?; 

inverting  the  order  of  integration  and  taking  account  of  (6.11)  and  (6.14),  the  re¬ 
sult  follows  immediately. 

As  an  application  of  (7.12)  let  us  consider 

L  M  x]=f(x)  =  g  (/S,  *  +  y)  ; 

if  we  remark  that  by  the  very  definition  of  the  Lebesgue  integral 

/>  /»  +  oa 

I  L[u\x]  dm  =  I  £d#i(Z), 

J  Q  •/  — 00 

we  get  immediately  the  condition  (6.7). 

8.  An  interesting  application  is  the  computation  of  the  statistical  average  of  the 
functional 

G[  u\x]  =  Y[Tta\xYP[Ttfa\xYV[Tta\xY\ 

according  to  (5.5),  (5.6)  and  (5.7),  this  average  gives  the  value  of  the  moment 

y  (%,  t) kp  (%,  t)  t)n. 

For  the  three  first  order  moments,  since  F,  P ,  V  are  linear  functionals,  we  have 
simply: 

(8.1)  y(x, t)  =  0, 

(8.2)  p(x,  t)  =  0  , 

(8.3)  v  (x,  t)  =  0  . 

In  order  to  compute  the  six  second  order  moments  we  must  introduce  a  function 
p(h )  defined  in  the  following  way: 

(8.4)  P(h,p)=YiJ^  g(P,  7)  g(P,  7  +  h)  dy , 
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(8.5) 
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p(h)  =  f  p  ( h ,  0)  dv. 

J  B 

We  shall  call  p(h,  0)  the  time  correlation  function  and  p(h)  the  statistical  correla¬ 
tion  function. 

It  is  quite  easy  to  prove  that: 

(8.6)  p  (h)  is  periodic  ,  p  (h  +  21)  =  p  ( h )  ; 

(8.7)  p  (h)  is  continuous  ; 

(8.8)  .  flp(h)dh  =  0; 

•'0 


(8.9)  p  (—  h)  =  p  ( h ) ; 

(8.10)  |p(A)|Sp(0); 

(8.11)  p  (h)  is  positive  definite6 ; 

+  °°  L 

(8.12)  p(A)  =  £c„cos»?f , 

1  1 


the  series  being  uniformly  convergent  and 

=  K(0) 2 +M« 

where 

a"(/3)  g  7)  cos  n  T  dy  ’ 

&»(/3)  =\f0  g(0>  7)  sin  n  dy. 

To  every  function  g  satisfying  (6.10)  to  (6.14) — that  is  to  every  convenient  choice 
of  an  invariant  measure  m  in  the  phase  space  12 — there  corresponds  a  statistical 
correlation  function  p(h)  with  the  properties  (8.6)  to  (8.12).  By  very  easy  com¬ 
putations  we  can  now  express  the  six  second  order  moments  as  follows: 

(8.13)  y(x,t)2  =  J  J  p(s2— Si)  dsidsz,  O^x^l, 

(8.14)  J{^1T2=2  [p(0)  +p(2*)], 


(8.15)  v  (x,  t)2  =  2  [  p  (0)  -p(2*)], 

(8.16)  y(x,  t)p(x,  t)  =\f  [p(s  +  x)  +  p(s  -  x)  ]  ds  , 

8  Positive  definite  means  that  for  every  set  of  values  hi,  hn  and  for  every  set  of  complex 
numbers  X\,  .  .  .  ,  Xn,  the  hermitian  form 

2  XiXtpihj  -  hk) 


is  nonnegative. 
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(8.17)  y  (x,  t)  v  (x,  t)  =  0, 


(8. 1 8)  p  (x,  t)  v  (x,  t)  —  0  . 

As  an  application  of  these  formulas  let  us  consider  the  functional  sometimes 
called  the  total  energy  of  the  vibrating  string: 

E\Tt «]  =  -2  flf{t  +  x)f(t-x)dx 
‘/0 

=  i  ty--p(x,  t)‘)dx; 

l  J  0 

we  have,  from  (8.14),  (8.15)  and  (8.8), 

(8.19)  E  [7\«]  =  0  . 

9.  As  we  have  seen,  the  time  average  and  the  statistical  average  of  a  nonlinear 
functional  are  in  general  different  in  the  statistical  mechanics  of  the  vibrating 
string. 

The  only  case  in  which  we  have  metric  transitivity  for  the  transformation  of  the 
phase  space  12  into  itself  corresponds  to  the  case  that  g(i 3,  y)  does  not  depend  on  /3: 

(9.1)  g(0,y)  =  g  (7)  • 

The  meaning  of  this  very  particular  choice  of  g  is  quite  simple:  in  this  case  our 
subset  12*  is  identical  with  a  given  trajectory  I\.  In  other  terms,  for  the  corre¬ 
sponding  measure  m  in  the  phase  space  12,  we  have 

tn  (TO  =  1  , 

and  for  the  set  12x  of  all  the  other  trajectories 

m  (12i)  =  0  . 

In  the  language  of  probability  we  assume  that  all  the  states  of  the  vibrating 
string  not  represented  by  points  w  belonging  to  the  particular  trajectory  Tx  have  not 
the  slightest  chance  of  being  observed. 

If  (9.1)  is  true,  for  every  functional  F[co]  €  L  we  can  put 

F  [«]  =  4>  (7)  ; 

then  (7.4)  becomes 

F[u]  =Ylf0  dy  ’ 
which,  according  to  (7.9),  gives  obviously 

(9.2)  Flw]  =MtF[ «]. 

10.  As  a  final  remark,  we  should  like  to  say:  if  someday,  starting  from  the  equa¬ 
tions  of  a  perfect  or  a  viscous  fluid,  the  moments  of  the  functions  involved  could 
be  computed  by  the  same  method  we  have  applied  to  obtain  (8.1)  to  (8.3)  and 
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(8.13)  to  (8.18),  then  it  would  be  all  right  to  speak  of  a  statistical  mechanics  of 
turbulence. 

What  has  been  done  so  far,  fruitful  and  important  as  it  has  been  for  the  develop¬ 
ment  of  fluid  mechanics,  is  a  quite  different  thing.  Since  no  invariant  measure  has 
been  introduced  in  the  phase  space  representing  the  states  of  a  fluid,  we  are  always 
concerned  with  more  or  less  empirical  definitions  of  some  kind  of  time  average; 
often  the  only  connection  of  these  averages  with  the  theoretical  equations  of 
fluid  mechanics  is  through  dimensional  analysis;  sometimes  a  really  powerful 
physical  intuition  seems  to  have  illuminated  the  whole  field  and  one  can  not  help 
admiring  the  advances  which  have  been  made.  But  from  the  theoretical  point  of 
view,  we  dare  say  that  only  the  first  steps  have  been  taken:  there  is  still  a  long  way 
to  a  statistical  mechanics  of  turbulence. 
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